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An efficient quantum algorithm is proposed to solve in polynomial time 
the parity problem, one of the hardest problems both in conventional quan- 
tum computation and in classical computation, on NMR quantum computers. 
\Q ', It is based on the quantum parallelism principle in a quantum ensemble, the 

selective decoherence manipulation, and the NMR phase-sensitive measure- 
QQ ' ment. The quantum circuit for the quantum algorithm is designed explicitly. 

o ■ 

1. Introduction 

f-| , It is well known that the parity problem is one of the hardest problems 

CL|| [1, 2]. These hard problems [1-4] including the parity problem, majority 

problem, and some iteration problem, etc., can not be solved efficiently both 

■ on the classical and quantum computers in polynomial time so far. It has 
^ • been shown that quantum computers can achieve no more than a polyno- 
mial speedup over the classical computers in solving these hard problems, 
implying that a quantum computer cannot outperform a classical computer 

^ I in solving these problems [1, 2]. On the other hand, since nuclear magnetic 

■ resonance (NMR) technique [5, 6] has been proposed to implement in ex- 
periment quantum computation [7, 8], there are some difficulties, limits, and 
questions in NMR quantum computation [9, 10], although a great and rapid 
progress has been achieved and a great flow of the work has been reported 
in past several years. Due to too low spin polarization in a spin ensem- 
ble at room temperature [9] it can not be certain up to now from quantum 
mechanical principle whether NMR quantum computation is a real quantum 
computation and can be more powerful than the classical computation. Some 
peoples [11-13] even claimed from the view point of the quantum entangle- 
ment theory that NMR quantum computers are not in fact real quantum 
computers, and suspect whether or not NMR quantum computation based 
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on the effective pure states [7, 8] is as powerful as the quantum computation 
based on the pure quantum states also owing to the two main drawbacks, 
that is, an exponent reduction of NMR signal intensity as the qubit num- 
ber and a high cost for the preparation of the initial effective pure state on 
NMR quantum computation. However, most recently an efficient quantum 
algorithm running on the NMR quantum computers is proposed to solve in 
polynomial time the unsorted quantum search problem [14] . It is well known 
that the conventional quantum search algorithm, i.e., the Grover algorithm, 
based on the pure quantum states [15] or the effective pure states [16] cannot 
solve the search problem in polynomial time, although it is a quadric speedup 
algorithm over the classical search algorithms. Therefore, the new quantum 
search algorithm shows that quantum algorithms on the NMR quantum com- 
puters are not only as powerful as on the conventional quantum computers 
based on the pure quantum states but also some quantum computations run- 
ning on the NMR quantum computer can be more powerful than running on 
the conventional one. This algorithm is based on the quantum parallelism 
principle in quantum ensembles, the selective decoherence manipulation, and 
NMR phase sensitive measurement [14]. In this paper I further show that 
there is an efficient quantum algorithm to solve the parity problem in polyno- 
mial time on the NMR quantum computers. It is still based on the quantum 
parallelism principle in quantum ensembles, the selective decoherence ma- 
nipulation, and NMR phase sensitive measurement. The present quantum 
algorithm to efficiently solve the parity problem in polynomial time supports 
powerfully the belief that the NP hard problems can be efficiently solved in 
polynomial time on quantum computers and quantum computers can out- 
perform eventually the classical computers. 

2. Quantum parallelism in a quantum ensemble 

The quantum parallelism is one of the most basic principles in quantum 
computation [17]. In principle, the quantum parallelism offers the possi- 
bility that quantum computation can be more powerful than the classical 
computation, although it is usually difficult to extract the desired computa- 
tional output from a quantum system in quantum computation directly by 
the quantum mechanical measurement [18]. The quantum parallelism may 
be understood simply through the process that superposition of a quantum 
system is acted on by a unitary transformation of quantum computation. 
Since a quantum system can be at any superposition which can be a linear 
sum of the usual quantum computational bases, while number and functional 
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operation in mathematics usually may be represented by the computational 
bases and the unitary transformation, respectively, in quantum computation, 
then performing quantum mechanically the unitary transformation of the 
functional operation on the superposition means really that the functional 
operation is performed in a parallel form with all input numbers represented 
by those computational bases of the superposition. Because quantum compu- 
tation is a reversible computational process [19] an arbitrary quantum state 
of the quantum system is equivalent unitarily to a computational basis or 
any other quantum state of the system, that is, any pair of quantum states 
of the system can be unitarily converted completely each other in quantum 
computation [20] when there are not quantum measurement and decoher- 
ence process during the quantum computation. Now, if the superposition is 
expressed as a unitary transformation acting on some usual computational 
basis, the above functional operation process really can start at the compu- 
tational basis instead of the superposition. However, this does not mean that 
there is not the quantum parallelism in the computational process even when 
the initial input state is not superposition. This is because superposition of 
the quantum system is usually created during the unitary transformation of 
the functional operation. Therefore, in fact the quantum parallelism is also 
characterized by unitary transformation besides superposition. Due to this 
fact sometime the initial input state is included in a quantum algorithm [21]. 
However, it is key important to separate the initial input state represent- 
ing in mathematics number and the unitary transformation representing the 
functional operation mathematically when the quantum parallelism principle 
is extended to a quantum ensemble and can be able to play an important 
role in quantum computation. In general, a molecule spin system and its spin 
ensemble composed of the molecule spin system obey the unitary dynamics 
governed by the same spin Hamiltonian in nuclear magnetic resonance spec- 
troscopy [5, 6]. This point is very helpful for quantum computation making a 
transition from the quantum spin system version to the spin ensemble version. 
It has been shown that any unitary transformation of quantum computation 
can be constructed with the spin Hamiltonian of a spin system [22]. Then 
the unitary transformation is still available and takes the same form when 
quantum computation is carried out in the spin ensemble with the same spin 
Hamiltonian of the molecule spin system. The significant difference between 
the spin system and its ensemble is that there are mixed states in addition to 
the pure quantum states in the spin ensemble. However, just due to the fact 
that a spin ensemble can be in a mixed state, the spin ensemble can provide 
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the quantum parallelism principle with a much larger space for it to bring 
into play in quantum computation when the unitary transformation and the 
initial input state are independent on each other in quantum computation. 
This is because the initial input state of the quantum computation may be 
any mixed state in addition to the pure quantum state in a spin ensemble. 

In a quantum system with n qubits there are N — 2^ usual computa- 
tional bases or N any other eigenbases in the A^— dimensional Hilbert space 
of the quantum system, and any superposition of the quantum system can be 
expressed as a linear sum of the N eigenvectors in the vector representation. 
In the matrix representation any superposition of the quantum system is still 
expressed as a linear sum of N linearly independent quantum states which 
are now represented by matrices, but this matrix representation just gives 
the same information of the quantum system as the vector representation 
provides. For example, the effective pure state in a spin ensemble [7, 8] pro- 
vides a faithful matrix representation for the pure quantum state of the spin 
system in the spin ensemble, but it does not give more information than the 
pure quantum state in quantum computation. However, any superposition 
of a spin ensemble, that is, any mixed state or any density operator of a spin 
ensemble must be expanded as a sum of any 4"^ linearly independent ma- 
trix bases in the 4^*— dimensional Liouville space of the spin ensemble [5]. A 
mixed state is not always able to be unitarily transferred completely into an- 
other mixed state in a spin ensemble, and this is quite different from the case 
in a quantum system in quantum computation. If any unitary transforma- 
tion and the initial input state are independent on each other in a quantum 
algorithm whose quantum circuit is constructed with the spin Hamiltonian 
of a spin system or its spin ensemble, any mixed state of the spin ensemble 
can be used as the initial input state of the quantum algorithm run on the 
spin ensemble. Then, the initial input state of the quantum algorithm can be 
selected from the larger 4"— dimensional Liouville space of the spin ensemble 
instead of the smaller 2"— dimensional Hilbert space of the spin system. This 
could be helpful to design efficient quantum algorithms to solve some hard 
problems on a spin ensemble. 

How to solve mathematical problems in a spin ensemble like in a quantum 
system? The well-known scheme is based on the effective pure state [7, 8], 
but according to the scheme the mixed states in a spin ensemble really do 
not play an important role in the quantum computation. A possible scheme 
proposed here is that the complete information of the problems to be solved 
first is contained or hidden in unitary transformation of quantum compu- 
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tation or its corresponding effective Hamiltonian and / or the mixed states 
in the spin ensemble, and the simpler and better one is that only the uni- 
tary transformation contains the complete information of the problem to be 
solved. By choosing the suitable initial mixed states one runs the quantum 
algorithm on the spin ensemble, and then extracts the desired information 
of the problem by NMR phase-sensitive measurement technique [5, 6]. The 
latter scheme has been used to find the efficient quantum search algorithm 
on NMR quantum computers [14]. 

3. Selective decoherence manipulation 

All the coherent and the noncoherent components (the longitudinal mag- 
netization and spin order components in nonequilibrium state) in a spin 
ensemble decay irreversibly as time and at the same time the spin ensemble 
returns to the thermal equilibrium state irreversibly, but their decay rates are 
quite different [5, 6]. For example, in a spin ensemble the longitudinal relax- 
ation time is always longer than the transverse relaxation time, and for some 
macromolecules the longitudinal relaxation time can be even hundred and 
thousand times longer than the transverse relaxation time. The important 
fact is that the decay rates for some coherent and noncoherent components 
in a spin ensemble can be manipulated and controlled externally. For exam- 
ple, the nonzero-order multiple-quantum coherences including the transverse 
magnetizations are much sensitive to the inhomogeneous external magnetic 
field but the longitudinal magnetizations are not [5, 6]. Then one can make 
use of the inhomogeneous external magnetic field, for example, a gradient 
magnetic field [6], to control the decay of these multiple-quantum coher- 
ences, so that at a short time these coherent components are cancelled but 
only the zero-quantum coherences and the longitudinal magnetizations and 
spin order components almost keep unchanged. Finally, the residual zero- 
quantum coherences can be cancelled by the zero-quantum dephasing pulse 
[24] or z-filter [5, 6, 25] so that only the longitudinal magnetization and spin 
order components are retained selectively in the spin ensemble. Such selec- 
tive decoherence manipulation has been used extensively in nuclear magnetic 
resonance experiments. In quantum computation the quantum measurement 
is also a selective decoherence manipulation. The quantum state measured is 
retained but any other states decay to zero quickly in a quantum system when 
performing quantum measurement. The importance is how to exploit in a 
positive manner the selective decoherence manipulation to help to find new 
polynomial-time quantum algorithms in quantum computation. In quantum 
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computation the final goal is to extract the desired information of the prob- 
lem to be solved by a quantum algorithm from the quantum system, but one 
need not obtain the complete information of the quantum system. There- 
fore, one may use properly the selective decoherence manipulation to cancel 
the undesired information but obtain valuable information which is suffi- 
cient to solve the problem. The spin ensemble provides a very suitable place 
to perform the selective decoherence manipulation for quantum computa- 
tion since nuclear magnetic resonance technique is well-developed at present. 
The massive quantum parallelism makes it possible for quantum computation 
outperforming classical one, but the desired computational output could not 
be obtained easily due to the limit of the quantum measurement principle 
in a quantum system. However, since the selective decoherence manipula- 
tion in a spin ensemble can be performed more conveniently, the undesired 
coherent and noncoherent components created in quantum computation in a 
spin ensemble could be cancelled and as a consequence, the desired compu- 
tational results could be detected conveniently by the NMR phase-sensitive 
measurement technique [14]. Here, the selective decoherence manipulation is 
a positive effect in quantum computation. Though the decoherence effect is 
usually a negative effect on quantum computation and need to be overcome, 
the selective decoherence manipulation could be useful to help finding new 
efficient quantum algorithms. For example, most recently an efficient quan- 
tum search algorithm [14] based on the selective decoherence manipulation 
is proposed on the NMR quantum computers. 

4. The quantum algorithm to solve efficiently the parity prob- 
lem 

The parity problem can be outlined briefly below. Given a function f{x) 
which is defined on the integers x: < x < N — 1 [N = 2"). For any x 
{0 < X < N — 1) the function f{x) can only take either or —1. The parity 
of the function f{x) can be defined as [1] 

N-l 

Parif) = n fix) (1) 

x=0 

where f{x) = ±1, a; = 0, 1, N — 1. The parity problem of the function f{x) 
is how to determine the parity function par{f) . Obviously, the parity function 

par(f) can only take either +1 or -1. To determine the parity problem in 
classical computation one needs to call N times the function f{x) so that f{x) 
can be known for each x value and hence the parity function par{f) can be 
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determined from Eq.(l). In quantum computation evaluation of the function 
f{x) over all the x values < a; < A^ — 1 can be performed simultaneously by 
only one call of the functional operation of f{x) on the superposition of the 
quantum system due to the quantum parallelism, but the final result, that 
is, all values of the function f{x), can not be extracted simultaneously due to 
the limit of quantum mechanical measurement. As a result, the cost to solve 
the parity problem on a quantum computer is not really much lower than on 
a classical computer, implying that quantum computers can not outperform 
classical computers in determining the problem. However, the parity problem 
is a global-type or collective-type problem but not an individual one of the 
function f{x). Perhaps people could not evaluate each value of the function 
f(x) over all x values but can determine the parity function (1). 

First of all, the unitary transformation corresponding evaluation of the 
function f{x) is introduced. This unitary operation is the oracle unitary 
operation in the quantum algorithm to solve the parity problem. Assume 
that the evaluation of the function f{x) can be expressed by performing 
the oracle unitary operation Uf on any quantum state |a;)|'S') in a quantum 
system, 

Uf\x)\S) = f{x)\x)\S) = expi-t7rg{x))\x)\S) (2) 
where the auxiliary quantum state 15*) must be chosen properly [1, 14, 18] 
so that equation (2) holds, and the function g{x) is defined by f{x) — 
exp{—i7rg{x)), that is, g{x) = if f{x) — 1; g{x) = 1 if f{x) = —1. When 
the oracle unitary operation Uf acts once on any superposition of the quan- 
tum system all the values of the function f{x) over all a;: < a; < — 1 are 
evaluated simultaneously, 

N-l N~l 

Uf E a.,\x)\S) = E aJ{x)\x)\S). (3) 

x=0 x=0 

However, these function values are not simultaneously observable in the quan- 
tum system due to the limit of the quantum mechanical measurement. It has 
been shown that to determine certainly the parity problem on a quantum 
computer one needs to execute around N/2 times the oracle unitary opera- 
tion Uf [1, 2]. Equation (3) shows that the oracle unitary operation Uf is 
independent upon any quantum state \x) of the work qubits / (hereafter the 
symbol / denotes the work qubits) in the quantum system but related to the 
auxiliary quantum state \S) (hereafter the symbol S denotes the auxiliary 
qubits). 

In order to construct the explicit equivalent form for the oracle unitary 
operation U f one may use the selective phase-shift operations CgiO) to express 
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C^{7ig{x))\r)\S) = exp{-ing{x)6rx)\r)\S) = { 



the oracle unitary operation Uj. The selective phase-shift operation Cs{0) is 
a diagonal unitary operator of the longitudinal magnetization and spin order 
operator subspace and defined as [14, 22] 

Csi9) = expi-iODs) (4) 
where the diagonal operator Dg is defined as Dg — diag{0, ...,0,1,0, ...,0), 
that is, {Ds)ss = 1 for index s and {Ds)rr — for any other index r ^ s. 
Note that when the selective phase-shift operation Cx{0) acts on any quantum 
state |r)|S') a phase shift exp(— i^^) is created only when the index r — x, 
otherwise the quantum state |r)|5') keeps unchanged, 

\r)\S) iir^x 
f{x)\x)\S) ifr = x' 
and also note that any pair of the selective phase-shift operations commute 
each other. It can be seen from Eqs.(2) and (3) that the oracle unitary op- 
eration Uf is really equivalent to the product of all selective phase-shift 
operations Cx{T^g{x)), x — 0,1, N — 1, 

N-l 

Uf=Yl aA^g{x)). (5) 

x=0 

In a more general case, in the present quantum algorithm the oracle unitary 
operation [// of Eq.(5) can be replaced with the more general oracle unitary 
operation Uo{0) which is the product of the selective phase-shift operations 
Cx{Og{x)), X = 0,1, N — 1, with any phase angle 9 by choosing properly 
the auxiliary quantum state \S) = |0)|1) [14] 

N-l 

Uo{e)^UfV{e)Uf= n c^{eg{x)). (6) 

Note that the function g{x) = 0, 1 for f{x) = 1, —1, respectively, as shown 
in Eq.(2). The parity function par{f) of Eq.(l) can be expressed as 

A^-l N-l N-l 

par{f) = n f{x) = n exp(-i7r^(a;)) = exp(-i7r E g{x)). (7) 

x=0 x=0 x=0 

Define the integer phase parameter G as 

G='e9(x). (8) 

x=0 

Obviously, par{f) = +1 if the integer phase parameter G is an even number 
and par{f) = —1 if G is an odd number. Therefore, the parity problem is 
really determined completely by the parity of the integer phase parameter G. 
Obviously, the integer phase parameter G is a global or collective quantity. 
Its parity may be determined without knowing each individual value of the 
function f{x) and hence the parity problem of the function f{x) could be 
solved efficiently without knowing all the functional values f{x) in advance. 
In the above the oracle unitary operations Uf oi the functional operation 
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f{x) and Uo{0) are introduced in a quantum system. Now to extend these or- 
acle unitary operations to the ensemble composed of the quantum system the 
auxiliary quantum states \S) in the quantum system must be kept unchanged 
or be replaced by their corresponding equivalent effective pure states [23] in 
the quantum ensemble because the oracle unitary operations are related to 
the auxiliary quantum states \S), as can be seen in Eqs.(2), (3) and (6), but 
the quantum state of the work qubits / can be replaced by any mixed 
state or density operator of the ensemble [14]. 

To describe conveniently the evolution process of the quantum spin sys- 
tem or its spin ensemble under the action of the oracle unitary operation 
Uo{0) (6) one needs to exploit the following general unitary transformation 
when the selective phase-shift operation Cs{6) acts on any density operator 
Pj{0) of the work qubits I of the spin ensemble [26]: 

Cs{e)pi{o)Cs{e)-^ = pM - (i - cos^)[p,(o), d^u 

+i sin e[pj{0), Ds] + [{1- cos 9y + sin^ e]DsPi{0)Ds (9) 
where the commutation [pj(0), Ds]+ — pj{0)Ds + DsPj{0) and [pj{0),Ds] — 
Pj{0)Ds — DsPj{0). With the help of the unitary transformation (9) and the 
explicit expression (6) of the oracle unitary operation Uo{0) one can obtain 
easily the following unitary transformation when the oracle unitary operation 
Uo{9) acts on any density operator pj{0)\S){S\ of the ensemble 

uMPimoio)-' = pM 

N-l N-l 

-(l-cos^)[p,(0), E gis)D,]+ + isme[pj{0), E g{s)D,] 

s=0 s=0 
N-lN-l 

+ [(l-cos^)2+sin2^] E E g{s)g{t)D,pj{0)Dt. (10) 

s=0 t=0 

Here, without confusion the auxiliary pure quantum state or effective pure 
state 15*) (5*1 is omitted. Equation (10) is the basic unitary transformation 
to analyze the quantum algorithm to solve efficiently the party problem on 
an NMR quantum computer. The initial density operator p/(0) must be 
related to each quantum state |a;) {x = 0,1,.., N — 1), otherwise it keeps 
unchanged under the action of the selective phase-shift operation Cx{0), 
C:,;{9) p I {0)Cx{6)~^ = Pi{0)- and as a result, one may not determine certainly 
the parity function of Eq.(l). One of the proper initial density operators in 
a spin ensemble may be 

Pl{0)^f: Skhy (11) 
k=l 

where is the spin polarization parameter of spin k and the term pro- 
portional to the unity operator is omitted without losing generality. It has 
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proved that the density operator pj{0) of Eq.(ll) are related to all the quan- 
tum states \x) of the spin system in the spin ensemble [14]. Furthermore, it 
is convenient to employ the unit-number vector representation {a|} [14] to 
analyze the unitary transformation (10). In the unit-number vector repre- 
sentation the diagonal operator Dg can be expressed as [14] 

Ds = (|^i+af/i,)(8)...0(i£;fe+a|7fe,)0...(8)(i^„+a^7„,). (12) 
Assume that the work qubits / are prepared initially at the density operator 
of Eq.(ll) and the auxiliary qubit S is in the quantum state |>S')(S'| [23] in the 
spin ensemble. Then make the oracle unitary operation of Uo{0) on the spin 
ensemble. The spin ensemble will evolve according to the following unitary 
transformation obtained easily from Eq.(lO): 

n n 

Uo{0) ^khyUo{0)~ —J2 £khy 

k=l k=l 

N-1 n 

-{l-cosO) E gis) E (|^i + aUi.)0- 

s=0 k=l 

<S){lEk-i+al_Jk-iz) ®{£khy) ®{\Ek-\-i+al^ih+iz) (8) ■■■ 0(i-£^n+<4^) 

7V-1 n 

-Sine E g{s) E {\E^ + a{h,)®... 

s=0 k=l 

lS){lEk-i + al_Jk-iz) (^{ekalhx) <S){^Ek+i + afe+i4+u) ... 0(|-Bn + 

N-1 n 

+ [(l-cos^)2+sin2^] E 9{s)g{t) E [|(l+afa*)£;i + i(af+a*)7i,]0... 

t>s=0 k=l 

0[|(1 + 4_i4_i)Efc_i + i(a|_i + 4_i)4_u] 0[|efc(l - alaDhy] 
(8)[i(l + 4+i4+i)£;fc+i + |(a|+i + a*+i)7fc+ J ... 

(8)[i(l + «)En + U< + OU- (13) 
The mixed state of Eq.(13) of the spin ensemble after the action of the oracle 

unitary operation of Uo{9) is quite complicated, but it can be simplified by the 
selective decoherence manipulation. First, a hard 90 degree pulse along y-axis 
direction is applied to all the spins of the work qubits 7 of the spin ensemble 
and then a purge pulse [5, 6, 24, 25], that consists of the z-gradient field 
used to cancel all the nonzero multiple-quantum coherences including single- 
quantum coherences, and the zero-quantum dephasing pulse [24] used to 
cancel all the zero-quantum coherences, is applied to cancel all the undesired 
multiple-quantum coherences but only the longitudinal magnetization and 
spin order components keep unchanged. Then the final density operator is 
reduced to the form 

n N—l n 

Pf = {2/N) E Skhzi E g{s)al) = (2/iV) E SkPkhz (14) 

k=l s=0 k=l 



10 



where the phase angle 9 in the oracle unitary operation Uo{0) is taken as 90 
degree. The density operator of Eq.(14) determines the final observable 
NMR signal of any spin k {k = l,2,...,n) of the spin system in the spin 
ensemble, and the intensity of the NMR signal of the spin k is proportional 
to the sum Pk, which is an integer defined by 

7V-1 G 

Pk = {E 9{s)ai) =E 4 (15) 

s=0 t 

where the latter equality holds due to the fact that the function g{s) = 0, 1 
for any index s, and the second sum (the index t) runs over only those 
indexes s with g{s) = 1. It turns out easily that the integer can take any 
integer from —N/2 to N/2 since a\ = ±1 for any indexes k and t, and a 
half of all unity numbers {a^} take +1 for any given index k. It follows 
from Eq.(15) that the parity of the integers Pk {k = 1, 2, ...,n) is determined 
uniquely by that of the integer phase parameter G and in turn, the parity 
of the integer phase parameter G is also determined uniquely by that of 
the integer P^, that is, if the integer phase parameter G is an even number 
each integer P^ must also be an even number, and if G is an odd number 
each Pfc is an odd number too. Now, if every integer P^ is an odd number 
it follows from the density operator P/ of Eq.(14) that the NMR signal for 
any spin k {k — 1,2, ...,n) of the spin system in the ensemble always has 
a nonzero intensity. Obviously, if there is zero-intensity NMR signal of the 
density operator pj of Eq.(14) for any arbitrary spin k of the spin system in 
the ensemble it can be certain that all the integers {Pk} are even numbers 
(including zero) and hence the integer phase parameter G is also an even 
number. However, there is another possibility that the NMR signal of the 
density operator P/ of Eq.(14) for all the spins in the ensemble can have 
nonzero intensity even when all the integers {Pk} are nonzero even numbers. 
It is clear that the parity problem can be solved efficiently by the above 
process only when the NMR signal of the density operator (14) for one spin 
of the spin system at least in the ensemble has a zero intensity. If the NMR 
signal of the density operator of Eq.(14) for all spins in the ensemble have 
nonzero intensities the above process can not determine certainly the parity 
problem. To solve this problem one needs to modify the above quantum 
algorithm. First, a phase-shift operation Um{—'^/2) is introduced. This 
known phase-shift unitary operation is not an oracle unitary operation and 
is applied only on the work qubits / of the ensemble. It is defined by 

M 

C/M(-7r/2) =n CK-7r/2). (16) 
I 
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This phase-shift unitary operation is the product of M selective phase-shift 
operations Q(— 7r/2), while these selective phase-shift operations are given 
in advance and are required to meet some constraint that will be stated 
below. It turns out easily that the phase-shift unitary operator ?7m(— '7r/2) 
commutes with the oracle unitary operator Uo{9): [?7m(— 7r/2), t/o(^)] — 0. 
Now a modified quantum algorithm to solve the parity problem is given 
below. One first makes an oracle unitary operation of Uo{tt/2) on the initial 
density operator P/(0) of Eq.(ll) and then applies the known phase-shift 
unitary operation ?7m(— 7r/2) to the ensemble, then adds y-axis direction 90 
degree pulse 90° and a purge pulse just as the above process. It is easy to 
prove along the line of the above process that the final density operator now 
is given by 

Pf = (2/A^) E SkiPk - Mk)hz (17) 
^ k=i 

where the contribution of the part to the NMR signal of Eq.(17) comes 
from the phase-shift unitary operation Um{—'^/'^)- The integer is given 

by 

M 

Mk =E 4. (18) 

Note that unlike the integer Pk each integer Mk is given in advance and can 
take any integer from —N/2 to N/2. Now the intensity of the NMR signal of 
the spin k of the spin system in the ensemble is proportional to the subtrac- 
tion (Pf. — Mk) of the two integers Pk and Mk instead of the single integer Pk- 
The subtraction (Pk — Mk) can be positive, negative, or zero. Obviously, The 
subtraction [G — M) is certainly an even number if there is zero-intensity 
NMR signal of the density operator Pf oi Eq.(17) for any arbitrary spin k of 
the spin system in the ensemble. In the case the parity of the integer phase 
parameter G can be determined directly from the subtraction (G — M) since 
the integer M is known , that is, G is an even number if the integer M is 
an even number, otherwise G is an odd number. How to find the zero point 
of the subtraction (P^ — Mk) with the above quantum algorithms? This is a 
classical one-dimensional search problem to find a zero point in mathematics. 
At the first step the original quantum algorithm is run once, that is, mak- 
ing the oracle unitary operation of Uo{t^/2) on the initial density operator of 
Eq.(ll), then adding a hard 90° pulse on the work qubits / and a purge pulse 
on the ensemble, then applying a readout 90° pulse on the work qubits / and 
recording the NMR signal of work qubits / during decoupling the auxiliary 
qubits S [14] by NMR phase-sensitive measurement [5, 6]. If one can find 
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that there is a zero-intensity NMR signal of any spin k in the fast Fourier 
transform NMR spectroscopy of the recorded NMR signal, it can be certain 
that the integer phase parameter G is an even number and the parity func- 
tion par{f) = +1, otherwise one needs to run further the above modified 
quantum algorithm. Without losing generality, choose an arbitrary spin k of 
the spin system in the ensemble with Pk > (note that N/2 > Pk> —N/2), 
then a suitable phase-shift operation [7^ (— vr/2) with = N/4 is set up (if 
Pk < 0, Mk = —N/4). By running the modified quantum algorithm one can 
find either (P^ — Mk) 7^ or [Pk + Mk) = from the fast Fourier transform 
NMR spectroscopy of the recorded NMR signal of Eq.(17). If {Pk - Mk) = 
the parity problem is solved. If {Pk — Mk) > a new phase-shift operation 
Um{—t^/'^) with Mk = 3N/8 is constructed, and if {Pk — Mk) < another 
phase-shift operation Um{—t^/2) with Mk = N/8 is built up. Then running 
the modified quantum algorithm again one can solve the parity problem if 
{Pk — Mk) = 0, otherwise construct again a new phase-shift unitary operation 
t^M(— 7r/2), and then run the modified quantum algorithm again till one finds 
the zero point {Pk — Mk) = 0. Obviously, the above search process really uses 
the classical one-dimensional two partition method to search for a zero point. 
It needs n times at most repeat running the quantum algorithm to find the 
zero point {Pk — Mk) = (note that N — 2") by changing the number Mk- 
Therefore, one can solve efficiently the parity problem by repeating n times 
at most to run the above quantum algorithm. 

5. The explicit construction of the phase-shift unitary operation 



The phase-shift operations 7r/2) is a diagonal unitary operator and 

its effective Hamiltonian can be expanded in the longitudinal magnetization 
and spin order operator subspace [22]. It is the product of M selective phase- 
shift operations C;(— 7r/2), as shown in Eq.(16). Without losing generality, 
one can set conveniently the index k — 1 and M — \Mi\ in Eq.(18), then for 
all M indexes / in the sum of Eq.(18) the unity numbers {a{} take the same 
value, i.e., a{ = +1 (Mi > 0), or = -1 (Mi < 0). Then the phase-shift 
unitary operation Um{—t^/2) in the unity-number vector representation {a\} 
may be expressed as, respectively. 



UM{-Tr/2) 



Um{-7t/2) = exp[m/2 E (|^i + hz) 0(|^2 + 4^2.) 




(19a) 
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or 

M 

C/M(-7r/2) = exp[i7r/2 E (l^i - hz) 0(|^2 + 4^2z) 

0...(8)(iK + aUn.)]- (19&) 
Without losing generality only f/M(— '7r/2) of Eq.(19a) with a\ = +1 for 
all index / is constructed explicitly below. One simple construction for the 
phase-shift unitary operations Um{,—t^I'^) is given explicitly by 
[/i(-7r/2) = exp[i7r/2(iEi + h,) 0(1^2 + hz) - 

)0(|K + 4.)], 

U2{-n/2) = exp[?7r/2(|^i + h,) 0(i^2 + hz) ... 0(|Sn-i + 4-1.)], 
[/3(-7r/2) = exp[^7r/2(iEi + hz) 0(|^2 + hz) ... 0(|K-2 + 4-2.)] 

X exp[z7r/2(iSi + Ji,) 0(1^2 + 4.) ... 0(|^n-2 + 4-2.)] 

0(iE„_i - 4_i,) 0(iE„ - 4,)], 

C/2.(-7r/2) = exp[i7r/2(|£;i + hz) 0(1^^2 + 4.) - 0(i^n-fe + 4-fez)] 

4Ar/2-i(-vr/2) = exp[z7r/2(i^i + 4.)] 

X expH7r/2(|Ei + hz) 0(1^2 - hz) ... 0(|K - hz)], 

4iv/2(-7r/2) = exp[i7r/2(|£;i + 4.)] 
where the unity operator Ek is omitted without confusion in the direct 
product terms, for example, (^-Ei + 4.) ... 0(|4i-i + 4-iz) = (^-^i + 
hz) ... 0(|-£'n-i + 4-iz) -E-n. Obviously, besides the above construction 
there are a lot of other different constructions for the phase-shift opera- 
tions which can be used to search for the zero point of (4k — M^) in the 
present quantum algorithm. Because all these phase-shift unitary opera- 
tions are not the oracle unitary operation and are given in advance one can 
construct their polynomial-time quantum circuits with the suitable interac- 
tion Hamiltonian of the spin system. In all the phase-shift unitary opera- 
tions C/m(— 7r/2) the basic unitary operations are C/2fc(— 7r/2) = Cq~^{—tt/2) 
{k = 0, 1, n-1), where C^-^(^) is a nonselective phase-shift operation in 
the spin subsystem with n — k qubits of the spin system with n qubits and 
is defined as Q"'^(^) = exp(-i^L)J^"'=) (note that C^{e) = Co(^)) with the 
diagonal operator D^''' = {^Ei + hz) 0(|^2 + hz) ... 0(|^n-fc + h-kz), 
(A; = 0,l,...,n— l).Itis well known that the nonselective phase-shift opera- 
tion Cq~^{9) {9 — tt) is the basic unitary operation in the Grover algorithm 
[15] with 2'^~^ dimensional search space. It is easy to prove that any phase- 
shift operation 4m(— vr/2) (1 < Af < N/2 = 2"-'^) can be expressed as the 
product of a polynomial number of the nonselective phase-shift operations 
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Cq *^(— 7r/2) {k = 0, 1, ...,n — 1) and the single-qubit operations. As an ex- 
ample, the phase-shift operation Um{-tt/2) with M = 2*^ 2' 2° < 2*^+^ 
can be expressed exphcitly as 

Um{-7t/2) = C7r^(-7r/2)exp(-m/„_fc,)Co"-'(-7r/2)exp(-^7r4_;,) 
xCo(-7r/2) exp(i7r/n-ia;) exp(i7r/n-te) 
In general case, the number M can be expressed by the binary representation 
as 

M = 2^ + hk-i2^-^ + 6fc_22*^-2 + ... + hi2^ + 6o2° < 2^^^ 
where 6; = 0, 1; / = 0, 1, A; — 1. Then the phase-shift operation Um{—t^/2) 
can be generally written as 

Um{-t^/2) = Q-*^(-7r/2)exp(-i7r4_fcJCo"-*^+^(-6fe_i7r/2)exp(-i7r7„_fc+i,) 

xCo"-'=+2j-6fc_27r/2)exp(-z7r/„_,+2,x)Co"-'=+'(-6fc-37r/2)exp(-z7r4_,+3x)... 

xCo"^(-6i7r/2) exp(-i7r/„_ia;)Co(-6o7r/2) exp(i7r/„_ij,)... exp(i7rJ„_fc^) 
Therefore, the phase-shift operation Um{—t^/2) can be expressed as the prod- 
uct of {k -\- 1) basic phase-shift operations Cq^'^^\—h /2) {I = 0,1,..., A;; 
< A; < 71 — 1) and 2k single-qubit unitary operations ex.p{±iTrIn-k+ix) 
{I — 0, 1, A; — 1) at most. 

6. Discussion 

It has showed in previous sections that the parity problem, one of the 
hardest problems in quantum and classical computations, can be solved effi- 
ciently on NMR quantum computers in polynomial time. It is based on the 
quantum parallelism principle in a quantum ensemble, selective decoherence 
manipulation, and the NMR phase-sensitive measurement. Because the par- 
ity function is a collective quantity one need not find all the function values 
so as to determine parity of the function on NMR quantum computers. In the 
paper the collective quantity G (see Eq.(8)) that reflects faithfully the parity 
problem is contained in the oracle unitary operations Uf and is used to de- 
termine efficiently the parity of function. The main drawback of the present 
quantum algorithm to solve the parity problem is that the detectable NMR 
signal reduces exponentially as the qubit number. To determine certainly the 
parity function with the present quantum algorithm it is required that the 
detectable NMR signal-to-noise ratio of the spin ensemble on an NMR ma- 
chine must be as low as 2ek/N, as can be seen in Eq.(14) or (17). These are 
similar to the basic characteristic feature of the conventional NMR quantum 
computers based on the effective pure states [9]. To implement really the 
present quantum algorithm one may still pay attention to the NMR dynamic 
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range problem on an NMR machine [5, 6]. This problem may be solved by 
using the spin ensemble of hetcromiclcar spin system. The present quantum 
algorithm is polynomial-time one, while those quantum algorithms running 
on the conventional quantum computers built up with the quantum system 
[1, 2] or based on the effective pure states [7, 8] cannot solve efficiently the 
parity problem and cannot outperform the classical algorithms. As an exam- 
ple, suppose that an NMR quantum computer with 30 qubits is available now 
[27]. Then to solve the parity problem with 30 qubits the present quantum 
algorithm needs ~ 30 calls of the oracle unitary operation Uo{7i/2) or about 
60 calls of the oracle unitary operation Uf. However, it will need ~ 2^°/2 
(~ 5 X 10^) calls of the oracle unitary operation Uf on the conventional 
quantum computer and 2^° (~ 10^) calls of the oracle unitary operation Uf 
on a classical computer, respectively [1] . The difference of the call number of 
the oracle unitary operation is huge between the present quantum algorithm 
and the conventional or the classical one. 

In comparison with the former scheme [14] to solve efficiently the parity 
problem based on the spectral labelling method [23] the present quantum 
algorithm is better one and more useful in practice. This is because the 
present quantum algorithm can be run on the spin ensemble of the linear 
molecules and does not required that the auxiliary qubits S are coupled with 
all the spins of the work qubits /, and moreover, it can use the highly mix 
state of the spin ensemble instead of the effective pure state as the initial 
input state. 

One could enhance the NMR signal of the final density operators of 
Eqs.(14) and (17) by running a quantum algorithm consisting of a poly- 
nomial number of the oracle unitary operations and nonselective unitary 
operations instead of the present quantum algorithm. However, the decoher- 
ence effect in a spin ensemble limits the possible call number of the oracle 
unitary operation in the quantum algorithm. Thus, there is a maximum 
call number of the oracle unitary operation in the quantum algorithm or a 
maximum run time Tmax of the quantum algorithm in the spin ensemble due 
to the decoherence effect. The NMR signal will decrease simply but not be 
enhanced again as increasing the call number of the oracle unitary operation 
in the quantum algorithm if the call number is larger than the maximum one. 
This is because at that case running the quantum algorithm will take more 
time than Tmax, so that the decoherence effect becomes the dominant effect 
negatively upon the NMR signal. If the NMR signal is still too weak to be 
detected on an NMR machine when the quantum algorithm is run with the 
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maximum number calls of the oracle unitary operation, then other NMR sig- 
nal enhancement methods such as the polarization transfer techniques, etc., 
[5, 6, 27-29] could be used to enhance the NMR signal and reach the same 
goal. Thus, here such NMR signal enhancement methods could be an effec- 
tive scheme to overcome the decoherence effect on the quantum computation. 
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